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Abstract 

A Kalman Filter alignment algorithm has been applied to cosmic -ray data. We discuss the alignment 
algorithm and an experiment-independent implementation including outlier rejection and treatment of 
weakly determined parameters. Using this implementation, the algorithm has been applied to data 
recorded with one CMS silicon tracker endcap. Results are compared to both photogrammetry mea- 
surements and data obtained from a dedicated hardware alignment system, and good agreement is 
observed. 
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1 Introduction 



Todays tracking detectors in particle physics experiments consist of several hundreds up to ten thousands of in- 
dependent detector elements, which allow to measure charged particle trajectories with a single-point resolution 
of typically 10-50 |J.m [Lil, J|. This resolution is significantly better than the placement accuracy as achieved 
during construction, which typically is an order of magnitude larger The established method to determine the 
true detector element position and orientation is to use measurements from particles traversing the detector in-situ, 
minimizing the residuals of an appropriate track model. Several algorithms have been proposed and exposed to 
data iHllSl. In this note we discuss the implementation and application of an alignment algorithm based on the 
Kalman Filter f6l to tracks from cosmic muons recorded by an integration setup of one CMS Tracker Endcap [7|, 
detailing the algorithm and summarizing results from reference [8J. The algorithmic implementation is done in a 
portable, experiment-independent manner and easily allows application in other experiments as well. 



2 Kalman Filter alignment algorithm 

The Kalman Filter algorithm has been implemented mainly according to the formulae described in reference ||6l 
with some modifications. 

The vector rh of detector measurements on a particle track is described by a function / that depends both on the 
track parameters p and the aUgnment parameters a, 

m = f{p, a) + e, (1) 

where the measurement errors are described by e, which has a known covariance matrix C. If the function / is not 
linear in the parameters {p, a), it is linearized at a starting point {pq, ciq): 

m = f{po,ao) + H{p^po) + D{a-do) + e + 0{{p~po)^,{a-ao)^) (2) 
~ fiPo, ao) - Hpo - DSq +Hp + Da (3) 



where the Jacobians H and D are 

df df 
H^^{po,ao), D=—{pn,ao). (4) 

Typically, for the expansion point ao the design geometry, knowledge from the assembly, or a previous alignment 
is chosen, and po are the track parameters obtained with this geometry. 

The goal of the Kalman Filter algorithm is to minimize the track residuals, i. e. to minimize the objective function 

fobj = {m- f{p, a))'^C{m ~ f{p, a)) (5) 

for the given track sample. This is achieved processing tracks in sequence and updating parameters and covariance 
matrix after each track. The resulting update equations for the alignment parameters a and their covariance matrix 
£;are L6,^ 

a = a + ED'^Wim- c - Dd) (6) 
E ^ E~ {ED^Y ■ W ■ ED^ (7) 

Here, W and V are auxiliary matrices, given by ||9l 

W = V-^ -V-'^H{H'^V-^H)-^H^V-^, V = C + DED^ . (8) 



Alignment parameters a and covariance matrix E are updated with each track. The algorithm needs a starting 
point for the parameters a and their covariance matrix E, which can be set to expectations e. g. from assembly 
tolerances. If in doubt, larger initial uncertainties are preferred since too small values can bias the alignment 
because the covariance matrix has decreasing eigenvalues. Parameters corresponding to global degrees of freedom 
can be fixed by assigning tiny prior uncertainties to specific alignment parameters. 
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In order to avoid bias from mis-measured tracks in the first steps, deterministic annealing with a configurable 
geometric schedule can be applied. This is done by increasing measurement uncertainties in Equation ([8]): 

V = aik)C + DED'', a(fc) = |^^ for 1 < fc < n, 

1^ 1 for K > n 

Here, k denotes the number of the current track, b is the annealing factor which is applied to the covariance matrix 
in the first step, and n is the number of the track from which on the full information from the track is kept. 

3 Implementation 

The software implementation of the algorithm has been done in an experiment-independent way. However, the 
input data are obtained from a specialized software framework tailored to the experiment, usually providing pattern 
recognition, track reconstruction as well as the description of the geometrical layout of the detector elements. 
Therefore, an interface between the experiment specific and the experiment independent software was designed. 
The basic choice was made that the experiment specific implementation has to provide all the information that is 
necessary to compute Equations (j6])-([8]) in a persistent matrix format, plus some additional information. 

3.1 Experiment specific implementation 

In the experiment specific implementation, the first additional information to be provided and stored is the number 
of parameters to be aligned. Parameters need to be identified uniquely by an index. This is necessary to save space 
and speed up the computation. The reason is that the Jacobian D in Eq. (|4]) has only few entries different from 
zero, since a single track only crosses few detector elements. 

A rough selection of tracks suitable for alignment has to be done, judging on the current alignment parameters. 
Then, for each track, the experiment specific software has to supply the measurement m with its covariance C, 
the constant c, and the Jacobians H and D. The Jacobians are evaluated according to the current knowledge at 
non-optimal parameters po and ao 

Additional experiment dependent information can be supplied, like run and event numbers and the number of the 
track in the current event. This can be especially useful if some tracks are rejected as outliers in the alignment 
procedure, such that one can have a look at the corresponding events within the experiment dependent software 
(e. g. event display). Also the value of the track fit, defined as 

— r'^ C f with r = {'m — c) (10) 

and the number of degrees of freedom fidof in the fit are stored. One has to note that in the case of an unaligned 
detector, does not follow a -distribution, but still provides some power to discriminate bad tracks. 

3.2 Implementation of the Kalman Filter alignment algorithm 

The Kalman Filter alignment algorithm is implemented as a C-n- program, which uses the ROOT ifTOl data analysis 
framework. The program reads the needed track information from the files which were created with the experiment 
specific program. 

The program initializes a vector of alignment parameters and the corresponding covariance matrix with config- 
urable starting values. One option is to pass over this information from the experiment specific part. However, 
pattern recognition and track reconstruction are very costly time-wise, and innut values to alignment only change 
little (e. g. due to non-linearities of f{p, a), which are neglected in Equation (Bb) when a different starting point is 
taken. By choosing the initialization to happen in the experiment independent part, avoiding a new reconstruction, 
computing time is saved. 

Then, processing track after track, the parameters a and the covariance matrix E are updated as specified in 
Equations (j6]l and (|7]l. Memory and computing time are saved in this step by reducing the alignment derivatives 
matrix D to non-zero columns, which correspond to the detector elements hit by the current track. In the same 
spirit, computing time is saved by using only the (known) non-zero elements of when computing the product 

^' A reference implementation in C++ |11| for the CMS experiment is available and can be used as a template for other 
experiments. 
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ED^ . The matrix _E_D^ consists of one line for each alignment parameter (not only those hit by the current track), 
and one column for each single measurement of the current track. The most computing time intensive matrix 
operation is the update of the matrix E, which is quadratic in the number of parameters and dominates the time 
consumption of the alignment procedure in case of many alignment parameters. 

During the update, the current alignment parameters and their uncertainties (extracted from the diagonal elements 
of the covariance matrix) are filled into histograms. After all tracks have been processed, these histograms to- 
gether with the final parameters, their uncertainties, and additional information on the alignment procedure like 
the number of hits for each detector element are stored. 



3.2.1 Outlier rejection 

To prevent bad measurements (e. g. from noise clusters) from having a large impact on the alignment results, outlier 
rejection is implemented. As mentioned above, the value (cf. Eq. ( 10 1) is computed for each track. The track is 
used for updating a and E when the probability P{x^ ,nj,of) exceeds a configurable threshold value. 

An additional outlier rejection is implemented by rejecting tracks that change the alignment parameters by more 
than a configurable factor beyond their uncertainties. Tracks showing this effect are assumed to be wrongly mea- 
sured. Considering the fact that the current alignment parameters combine the information of all previously pro- 
cessed tracks, one additional track should not have such a large impact on the results. 



3.2.2 Weakly determined parameters 

The objective function ^ is invariant under a global translation and rotation of all detector elements, i. e. a special 
linear transformation of the alignment parameters. Even when not present in the starting values a and E, such 
transformations can build up due to mis-measurements, round-off problems in the alignment procedure, or an 
incomplete track model. Large alignment parameters a have a negative impact on the alignment procedure due to 
larger distance to the linearization point. Therefore, these global parameters should be fixed. This can be achieved 
by assigning small uncertainties to some (linear combinations of) parameters for the initial matrix E, depending 
on the experimental setup, which corresponds to the definition of a geometry reference system. 

Furthermore, large eigenvalues can be present in E even after processing the last track. These large eigenvalues 
correspond to certain linear combinations of alignment parameters which geometrically represent a systematic 
distortion of the detector units position and orientation that is only weakly determined. This happens especially if 
the track sample consists of tracks with similar topology. Artificial distortions can bias physics observables like 
invariant masses, momentum scale etc. 

We consider two methods to suppress weak modes: The first method is to assign a small initial error to the weakly 
determined linear combination of parameters in E. The alternative is to align without changing the initial values 
of E, and later fit the amplitude of the weakly determined mode and subtract it from the parameters. 



4 Application to data 

The Kalman Filter alignment algorithm was applied to data taken from the integration setup of one endcap of the 
CMS experiment ll3j|2l tracker (TEC-n). The integration took place in 2006 at RWTH Aachen University. Apart 
from commissioning the system hardware, tracks from cosmic muons originating from air-showers were recorded 
and used for alignment. 

4.1 CMS Tracker Endcap 

The CMS tracker ^ is entirely based on silicon detector technology. It can be divided into five subsystems 
(Fig. [T]i: Pixel detector (PIXEL), Tracker Inner Barrel (TIB), Tracker Outer Barrel (TOB), Tracker Inner Disks 
(TID-H, TID-), and Tracker Endcaps (TEC-n, TEC-). 

CMS uses the following coordinate definitions: The y-axis points upwards, the x-axis points radially to the centre 
of the LHC lfT2l ring, and the z-axis points in direction of the beam line, completing a right-handed coordinate 
system. The azimuthal angle cj) is measured to the a;-axis in the x-y-plane, and hence describes rotations around 
the z-axis, and the polar angle 9 is measured to the z-axis. 

The Tracker Endcap TECh- covers the range 0.9 < ?; < 2.5 and consists of 3200 trapezoidal silicon-strip detectors. 
Both TECs consist of nine disks (Fig. |2| carrying 16 substructures called petals. Eight petals are mounted on the 
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Figure 1: View of the CMS tracker in the rz-plane fi\. Each Hne in the strip tracker represents a siHcon strip 
detector, whereas lines in the pixel detector represent ladders and petals on which the detectors are mounted in the 
barrel and endcaps, respectively. 



side facing the interaction point (front petals) and eight on the far side (back petals). By grouping the neighbouring 
front and back petals, the disks are divided into eight sectors numbered from 1 to 8. On the petals (Fig.[3]l, detectors 
are mounted with the strips in radial direction in up to seven rings. 




Figure 2: Left: Sketch of a tracker endcap O. Middle: Photo of a TEC disk. Right: Sector numbering scheme. 



4.2 Sector tests 

Petals in TEC+ were integrated and commissioned sector by sector The petals were mounted and subsequently 
attached to services (cooling, power, trigger, communication and data lines). Connections were tested with custom 
hard- and software. Noise data were recorded in order to spot potential flaws, allowing for repair or replacement. 
As a final test, signals from cosmic muons originating from ak-showers and traversing the detector were recorded 
and used for various studies, including alignment QO. 

Figure |4] shows the experimental setup. During cosmic muon data taking, TEC+ was in a vertical position (disk 
planes horizontal). Two square areas of size 80 cm x 80 cm, one below and one above the TEC+ sector under 
study, were covered with four AMS 1 13| scintillator panels 1 14|, equipped with two photo-multipliers on opposite 
sides. A coincidence signal from one of the upper scintillator panels and one of the lower panels was demanded 
to trigger a readout of the TECh- sector. A 10 cm thick lead shield was placed below TEC-n, but above the lower 
scintillators, in order to absorb low energy (< 250 MeV) muons and prevent them from triggering readout. 

Customized CMS software was used to read out the detector, detect signals, and reconstruct tracks. Table[T]shows 
the number of triggers and reconstructed tracks used for alignment for each sector in chronological order. 
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Figure 3: A TEC front petal (left) and a back petal (right) fl]. Only detector units on the side facing the interaction 
point can be seen (rings 1, 3, 5, and 7). Detector units on rings 2, 4, and 6 ai^e located on the back side of the petal. 



4.3 Results 

4.3.1 Track selection 

Tracks were reconstructed with bloated hit uncertainties in order to allow an efficient track reconstruction in the 
presence of misalignment and selected with loose criteria. To account for the overestimate, tracks with a x^- 
probability P{x^, n^of) < 0-5 with from Eq. 



10 1 were rejected. 



After this selection, 4000 to 6000 tracks were available for each of the sectors 1, 3, 5 and 7. Due to an optimized 
trigger configuration, which took place only after integration of the odd-numbered sectors, there were between 
35 000 and 60 000 tracks available for the even-numbered sectors 2, 4, 6 and 8. In total, the sector test data contain 
about 220000 tracks which were used in alignment studies. 

The starting geometry for the linearization of Eq. ([3]) was the design geometry, and the resulting measurements, 
predictions, derivatives, and covariance matrices were stored. 



4.3.2 Disk alignment 

Higher level structures like TEC petals and disks were aligned. Here, we describe the alignment of TEC disks, 
for which determined alignment parameters can be compared directly to measurements from survey and the Laser 
Ahgnment System (LAS). 

Each of the eight TECh- sectors has been tested separately. Therefore, no information can be gathered from the 
collected data about the relative position of the sectors with respect to each other TECh- disks were aligned by 
using all available tracks and assuming that detector elements were mounted on their nominal position on the petals, 
and petals on their nominal position on the disks. For each disk, corrections Ax, Ay, and Acj) were calculated. The 
corrections to the geometry have the same absolute value as the alignment parameters in Eq. (j6]l, but with inverted 
sign. 

The Kalman Filter alignment algorithm was configured to use zero as the initial value of the alignment parameters, 
corresponding to design geometry. The covariance matrix was initialised with large startup errors of 10 cm and 
10 rad for spatial and angular parameters, respectively, to avoid biasing the alignment results by the starting values. 
The reference system was defined by assigning very small initial uncertainties of 10^^ cm and 10^^ rad to the x, 
y, and (p parameters, respectively, for the disk closest to the interaction region. A standard annealing configuration 
with 6= 10 000 and n — 100 was used. Tracks which resulted in an update of the parameters larger than the 
parameter uncertainty were rejected as outliers. 



Table 1: Sector test data sets in chronological order. 



Run 


Sector 


Number of 


Number of 


number 




triggers 


tracks 


20944 - 20952 


5 


87400 


6961 


21136-21163 


3 


122058 


5821 


21238 -21269 


1 


88591 


5807 


21428 -21448 


7 


83772 


8502 


21512-21530 


2 


82139 


62650* 


21592-21617 


6 


87258 


68210* 


21666-21670 


4 


58248 


42260* 


21713 


8 


73227 


61159* 



* optimized trigger configuration 
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To verify convergence, several starting points for the parameters were tried and the order in which the data were 
processed was changed from chronological to random order. In both cases no significant changes in the resulting 
alignment parameters was observed. 

Due to the low number of aligned parameters, the time used for computing the alignment constants, roughly 300 s, 
was dominated by input/output operation of the input data. 

Due to the nature of the problem, the objective function only changes minimally when the hnear transformations 

Aa;(z) = Ax' (z) + + bx ■ z 
Ay{z) = Ay'{z) + ay + by ■ z 
A(j){z) = A(j)' (z) + + ■ z 

are applied to the parameters, where z is the coordinate shown in Figure [T] These transformations correspond 
to weakly determined deformations of the TEC+ structure. The computed alignment parameters were, after the 
alignment procedure ended, transformed such that 

Axi ■ Zi = 

i 

^ Aa;, = 

i 

and corresponding constraints for the Ayi and A(/)i were fulfilled. 

Figure |5] shows the corrections Aa; for disks 3 and 9 as a function of the number of processed tracks. It can be 
seen that the uncertainty on the correction decreases with the number of processed tracks. The uncertainty is larger 
for corrections belonging to alignables which are farther away from the disk that is used as the reference system 
(disk 1). 

Table|2]lists the obtained alignment corrections together with their errors. The result is also displayed graphically 
in Figure |6] The size of the corrections is of order 100 |J.m in x and y, and 100 i^rad in which corresponds 
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to an arc length of about 100 ixm at the outer disk circumference. The accuracy of the positional and rotational 
corrections, 2-10 \xm and 2-10 |j,rad, respectively, is very high compared to the typical size of the correction. 



Table 2: Alignment corrections obtained from disk alignment in x, y. 



Disk number 


Ax [ \xm] 


Ay [\mi] 


A(f> [larad] 


1 


24 ± 8 


-99 ± 8 


33.6 ±7.3 


2 


60 ±6 


-77 ± 6 


45.3 ± 5.6 


3 


-20 ±4 


17 ± 4 


51.2 ±4.1 


4 


-26 ± 3 


76 ± 3 


-53.6 ± 2.7 


5 


83 ± 2 


146 ±2 


-117.0± 1.7 


6 


-162 ±2 


42 ±2 


-111.4± 1.9 


7 


-67 ± 4 


21 ± 4 


1 15.8 ± 3.7 


8 


7±6 


-51 ± 7 


29.9 ± 6.2 


9 


100 ± 10 


-74 ± 10 


6.2 ± 9.2 
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Figure 5: Corrections Aa; for disks 3 (left) and 9 (right) as a function of the number of processed tracks. The width 
of the band represents the uncertainty on the correction. 
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Figure 6: Alignment corrections for disk alignment in x, y, 4>. Error bars are smaller than the symbol size. 



4.3.3 Comparison to Laser Alignment System results 

The TEC also comprises a hardware alignment system ||3|, using sixteen laser rays of wavelength A = 1075 nm 
propagating parallel to the beam line. All rays pass through the back petals, eight rays pass at a radial distance 
r4 = 564 mm and eight further rays at rg — 840 mm, the subscript indicating the ring number. The rays pass 
through an opening in the silicon detector's backside metallization. The laser rays get partially absorbed and 
produce signal on the detector elements. Using a straight line hypothesis, residuals can be deduced. These residuals 
were used to determine corrections to position and rotation of the nine TEC disks. Ax, Ay, and Acj) |15J. The 
corrections are listed in Table|3]and compared to the results deduced from the Kalman Filter in Figure]?] 

The RMS of the differences between Laser Alignment System and Kalman Filter corrections is 70 \xm in Ax, 
52 |j.m in Ay, and 56 [4rad in A(j>. Although overall good agreement is observed, some measurements deviate 
significantly. However, perfect agreement is not expected since the Laser Alignment System relies exclusively on 
144 measurements of detector elements in rings 4 and 6 of back petals only, whereas the track based alignment 
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Figure 7: Comparison of Laser Alignment System ifTSl and track-based alignment results. 



comprises data from more than 3000 modules on both front and back petals. In a further analysis it was found that 
the differences are compatible with the mounting precision of the petals. Especially the large deviation in Ax and 
Ay for the disk at z « 2050 mm could be traced back to an imperfect front petal mounting in sector eight. Here, 
one of the precision pins with which the petals are attached to the disks was visually confirmed to be bent into the 
direction indicated by the alignment corrections. 

Interpreting the result, the largest contribution to the misalignment is due to the positioning of the disks in the 
TEC-H, and the remaining difference is to a large extent due to the positioning of the petals on the disks. This is in 
agreement with the naive expectation that smaller structures can be assembled and mounted with higher precision 
than larger structures. 

4.3.4 Comparison to survey measurements 

A further validation of the alignment results was performed with survey data recorded by a survey team before 
integrating the petals and after their integration. Using photogrammetry, the displacements of four points at the 
outer circumference of each disk were measured with a precision of about 60 |J.m. Each measurement was done 
once with the TECh- being in horizontal and once with it being in vertical orientation. As before, corrections Aa; 
and Ay to the disk position as well as A(p to the disk rotation were estimated from these measurements. 

Table 3: Displacements Ax, Ay, A0 of TECh- disks determined with LAS residuals. The precision of the position 
corrections Aa; and Ay is 23 |a.m, and A0 is determined with an accuracy of 23 [a.rad. 



Disk number 


Ax [yim] 


Ay [yxm] 


A(j) [yxrad] 


1 


-64 


-6 


126 


2 


17 


-99 


30 


3 


7 


-33 


-47 


4 


43 


43 


-78 


5 


132 


89 


-116 


6 


-30 


116 


-38 


7 


-123 


-6 


60 


8 


-67 


-73 


34 


9 


86 


-31 


31 
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The measurements were transformed into the same reference system as used for the track-based ahgnment and are 
Hsted in Table |4] A graphical comparison of the track-based alignment results and survey data is given in Figure |8] 
The track-based alignment results agree with the survey measurements within (RMS) 63 [im in Ax, 52 [im in Ay, 
and 34 |xrad in A0 and show, within measurement precision, good agreement. 

Table 4: Displacements Ax, Ay, A0 of TEC+ disks determined with survey measurements. The precision of the 
position corrections Ace and Ay is 57 |J.m, and A^ is determined with an accuracy of 47 |a.rad. 



Disk number 


Ax [i^m] 


Ay [urn] 


A{j> [[orad] 


1 


30 


-63 


63 


2 


81 


-54 


19 


3 


-15 


-15 


57 


4 


18 


6 


-5 


5 


-31 


110 


-177 


6 


-73 


158 


-55 


7 


-165 


-21 


125 


8 


1 


-64 


20 


9 


154 


-56 


5 
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Figure 8: Comparison of survey and track-based alignment results. 



5 Summary 

A Kalman Filter alignment algorithm has been implemented in an experiment-independent way and applied to 
data. The obtained alignment corrections have been vaUdated with both a hardware alignment system and survey 
data. The difference (RMS) between the corrections estimated with the track-based alignment and the hardware 
alignment system is 70 and 52 [im in Aa; and Ay, respectively, and 56 [a.rad in A(f). The track-based alignment 
results agree with the survey measurements within (RMS) 63 [im in Ax, 52 [a.m in Ay, and 34 |j.rad in Acf). The 
difference betweeen laser alignment and survey data consistently is of a similar magnitude. Statistical uncertain- 
ties on the obtained alignment corrections are negligible when compared to the statistical uncertainties of laser 
alignment and survey data. 
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